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ABSTRACT 

The 'Bloated Stars Scenario' proposes that AGN broad hne emission originates in the 
winds or envelopes of bloated stars (BS). Alexander and Netzer (1994, 1996) estab- 
hshed that ~ 5x10^ BSs with dense, decelerating winds can reproduce the observed 
emission line spectrum and line profiles while avoiding rapid collisional destruction. 
Here, I investigate a third prediction of the model, related to the size of the broad line 
region, by deriving the emission line response to variations in the ionizing continuum 
('line reverberation') and comparing it to observations. The expected time lags, as well 
as the order of response of the various lines, strongly depend on the typical variability 
time scale of the ionizing continuum. The BS model studied here corresponds to a 
bright Syfert 1 galaxy (Syl) or a low luminosity QSO. I find that the BS model is 
consistent with the observed correlation between the Balmer lines time lags and the 
AGN luminosity, which at present is the only line reverberation information available 
for this luminosity class. The model displays also some of the anti-correlation between 
the time lags of the metal emission lines and their degree of ionization that has been 
observed in a few low-luminosity Syls. However, the bright Syl model results differs 
from the low-luminosity Syl data in that the H/? time lag is relatively shorter and the 
CivA1549 time lag relatively longer. 
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1 INTRODUCTION 

The observed properties of active galactic nuclei (AGN) lead 
to the conclusion that the broad line emission originates in 
numerous small, cold and dense gas concentrations, which 
are photoionized by the central continuum source. These 
objects are labeled "clouds", although their true nature re- 
mains unknown. A basic challenge for any broad line region 
(BLR) model is to specify the physical mechanism that pro- 
tects the clouds from rapid disintegration in the AGN's ex- 
treme environment, or else specify a source for their contin- 



ued replenishment. The bloated stars (BSs) model (Edward; 



198C|; [Mathews 1983t pcoville fc Norman 1988| ; [Penston 1988t 



Kazanas 1989) proposes that the lines are emitted from the 



winds or mass loss envelopes of giant stars. The star provides 
both the gravitational confinement and the mass reservoir 
for replacing the gas that evaporates from the envelope to 
the interstellar medium (ISM). The BS model is further mo- 
tivated by the lack of observational evidence for net radial 
motion in the BLR (e.g. Maoz et al. 1991, Wanders et al. 
1995). This is consistent with virialized motion in the grav- 
itational potential of the nucleus. 
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The BS model was recently studied by Alexan der fc 
Netzer ( |l994| ) (paper I) and Alexander & Netzer ( |l99e| ) 
(paper II), which used a detailed photoionization code to 
calculate the emission line spectrum and profiles of various 
BS wind models. The integrated emission line spectrum was 
calculated by combining the line emissivity of the BSs with 
theoretical models of the stellar distribution function. The 
model results were compared to the mean AGN line ratios 
and to estimates of the BLR size and line reddening. The 
mean observed Lya equivalent width was used to determine 
the number of BSs and their fraction in the stellar popula- 
tion and to estimate the collisional mass loss rate and its 
effect on the ISM electron scattering optical depth. Detailed 
models of the stellar velocity field were then used to derive 
the line profiles, which in turn constrained the BS distribu- 
tion within the stellar population. 

The main result of papers I and II is that there are BS 
wind models that can reproduce the emission line ratios and 
profiles to a fair degree with only a small fraction of BSs in 
the stellar population. However, these models do not resem- 
ble winds of normal supergiants. The photoionization calcu- 
lations show that the emission line spectrum is dominated 
by the conditions at the outer boundary of the line emitting 
zone of the wind. The successful BS models are those with 
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dense envelopes that have small density gradients. In this 
case, the wind boundary is set by tidal forces near the black 
hole and by the finite mass of the wind at larger radii. Only 
~ 5x10'' such BSs (< 1% of the BLR stellar population) are 
required for reproducing the BLR emission. As a result, the 
collision rate is reasonably small (< IMq/jx). On the other 
hand, lower-density wind models or those with steep density 
gradients are ruled out because they emit strong broad for- 
bidden lines, which are not observed in AGN. In addition, 
their low line emissivity requires orders of magnitude more 
BSs in the BLR. This leads to rapid coUisional destruction, 
very high mass loss rate and very short BS lifetimes. Even if 
the BSs are continuously created, the ISM electron scatter- 
ing optical depth is very large, contrary to what is observed. 

Paper II shows that in order to reproduce the typical 
shapes and widths of the observed emission line profiles it 
is necessary to assume that the BS fraction in the stellar 



succeeded in reproducing the line ratios and profiles of a 
typical AGN, and here I check whether it is also consistent 
with the limited available line reverberation data. Section ^ 
briefly summarises the main components of the BS model. 
Section |^ describes the line reverberation simulation proce- 
dure. The simulated results are presented in section ^ and 
discussed in section H. 



2 THE BS MODEL 

As described in detail in papers I and II, t he BS model is 
specifled by two main ingredients. The stellar distribution 
around the black hole and the density structure of the BS 
envelope. The BS model that is discussed in this paper is 
model C of paper II. 

The stellar distribution function that I use is based 



population falls off as 



The emerging picture is that on the numeric results of Murphy, Cohn & Durisen (1991) 



AGN harbour an exotic class of giant stars, whose existence 
is linked to the uniquely violent environment of AGN, and 
are consequently much rarer in the more usual stellar envi- 
ronments outside the inner nucleus. There are at least two 
astronomical precedences for such a situation: The existence 
of blue stragglers in the high density cores of stellar clusters, 
which is thought to be linked to stellar mergers or collisions 
(e.g. Bailyn 1995) and the recently discovered concentration 
of extremely rar e He I blue giants in the inner 1 pc of the 
Galactic Center (Krabbe et al. 1991). The envelopes of these 
blue giants reach densities as h igh as those assumed in the 
BS model ( [Najarro et al. 1996| ). 

The present work explores a third aspect of the BS 
model: BLR line reverberation, or echo mapping. The con- 
tinuum luminosity of AGN is observed to vary in time and 
in many cases, the continuum variability pattern is echoed 
after some time lag by the line fluxes. This is considered 
the most direct support for the photoionization hypothe- 
sis of BLR emission. The time lag between the continuum 
and line variations is interpreted as the difference between 
the light travel time along the direct line of sight and that 
along the indirect path from the continuum source, via the 
line emitting gas, to the observer. The continuum and line 
light curves are therefore related by the spatial emissivity 
distribution ('emissivity map') of the line, which in turn de- 
pends on the distribution of line emitting objects and their 
ionization structure. 

Extracting information on the BLR geometry from line 
reverberation data is severely limited by observational diffi- 
culties and the ensuing ambiguities in the statistical analysis 
(e.g. Maoz 1994). A reliable analysis requires good sampling 
over a long period, a condition which is seldom met. Only 
a few monitoring campaigns were carried out to date, each 
with a different duration and sampling frequency. In con- 
trast to the case of the line ratios and profiles, very little 
can be determined at this stage about the statistics of the 
time lags and their properties. However, some general trends 
are beginning to emerge: the BLR size seems to be positively 
correlated with the continuum luminosity, and high ioniza- 
tion metal lines see m to originate at s maller radii than the 
low ionization ones (Clavel et al. 1991). 

In this third part of the work I study the line reverber- 
ation properties of a BS model with a luminosity of a bright 
Seyfert 1 galaxy (Syl) or a low-luminosity QSO. This model 



(MOD), which follow the evolution of a spherically sym- 
metric multi-mass coeval stellar cluster in the presence of a 
central black hole. The black hole mass grows as it accretes 
mass from the stars, both directly, by tidal disruption, and 
indirectly, from mass loss in the course of stellar evolution 
and of inelastic stellar collisions. While the calculations deal 
in a self-consistent way with the stellar dynamics and the 
build-up rate of the central mass, the resulting AGN lumi- 
nosity is not well defined. This is because the calculations 
of the nucleus evolution are insensitive to the properties of 
the central mass, such as the ratio between the black hole 
mass and the infalling mass or to the black hole's angu- 
lar momentum, nor are they affected by the details of the 
accretion mechanism or by the rest-mass to luminosity con- 
version efficiency, rj. MCD make the simple assumption that 
the bolometric luminosity is I/boi(i) = riM{t)c? with r] = 0.1, 
where M is the total stellar mass loss rate from the stellar 
system. In this work I use MCD model 2B, at to = 10^ yr, 
when the black hole mass is Mbh = 1.9x lO*Af0, and the ion- 
izing luminosity -Lion = 3.6x10'*'' erg s~'. This luminosity 
is sub-Eddington and is determined by the total mass loss 
rate from the stellar system. This evolved system represents 
a luminous Syl or a low-luminosity QSO. Once the choice 
of the stellar model is made, the remaining free parameters 
are /as , the fraction of BSs in the normal stellar population 
and the inner and outer radii of the BLR, rin and rout. In 
paper II, it was shown that these parameters are constrained 
by the observed emission line profiles to values of roughly 
Tin = 0.001 pc, Tout = 0.25 pc, and /bs oc r~^, with the BSs 
comprising less than 1% of all stars. 

A single BS is modeled by a simplified, spherically 
symmetric structure with no claim to hydrodynamical self- 
consistency and without specifying the process that drives 
the wind. It has two components: a giant star of radius 
Ri, — lO'"^ cm and a mass of M* = O.SM©, emitting a spher- 
ically symmetric wind at a rate of Mw ~ Mq /yr. The 
wind, or envelope, extends to a radius J?w which contains an 
additional mass of up to Aiw = O.2M0. The wind density is 
parameterized as a power-law, with two free parameters 



N{R) = N^R/R^y 



(1) 



where A*'* is the hydrogen number density at the base of 
the envelope. The combined constraints of the observed line 
ratios and profiles flx these parameters to values of roughly 
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/3 = 3/2 wdN^ = 3.8x10" cm~^. The proximity of the BSs 
to the black hole makes tidal disruption the main physical 
mechanism that limits the size of the wind. For example, at 
30 Id, the BS wind radius (model C) is R„ ~ 2x10^^* cm. 

The BS model is fully specified by the choice of the AGN 
dynamical age, to, which fixes Mbh, by the BS structure and 
the ionizing continuum spectrum. However, tidally-limited 
BSs obey certain approximate scaling relations that link A^* 
and Lion in a way that allows Lion and A'^* to be jointly mod- 
ified without affecting the emissivity map (paper II). Since 
the emissivity map determines the line ratios and profiles, 
this scaling property means that the model luminosity can- 
not be currently fixed by either the BS model or the MCD 
AGN model. This introduces some uncertainty in matching 
the BS model with observed AGN that complicates the com- 
parison of the model time lags with the observed ones (sec- 
tion H). 




100 200 300 400 500 600 

Time (days) 



3 CALCULATIONS 

In order to calculate the light curves of the various emis- 
sion lines it is necessary to specify the line response to the 
local changes in the ionizing flux. A commonly assumed sim- 
plification is that the response is linear, in which case the 
relation between the line light curve, L{t), and the ionizing 
continuum light curve C{t) can be described by a convolu - 
tion with a transfer function ^' (Blandford fc McKee 1982) 



m = 



C{T)<l'{t~T)dT. 



(2) 



In the c ase of a spherically symmetric BLR the transfer func- 
tion is ([Maoz et al. 1991 ) 



r 0, 



■fr°n* w(r)c/2rdr, 
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< t < 2nn/c 

2nn/c< t <2rout/c 
2rout/c < t 



(3) 



The weight w{r) is the normalized line emissivity from a 
shell of BSs at radius r. 



/'"out 
r'^nBsir)Le{r)dr , (4) 

where Li{r) is the BS emissivity in line ^ at a distance r from 
the continuum source and hbs is the BS spatial density. w{r) 
is calculated by the photoionization code for the specific BS 
model. The assumption of a linear response implicitly in- 
cludes the approximation that the varying continuum does 
not change the emissivity map of the BLR, i.e. that \E' is 
unaffected by C{t). While this is not strictly true (the ion- 
ization parameter is proportional to the ionizing flux), it is 
probably an adequate approximation as long as the varia- 
tions are not large. Another way of interpreting '^{t) is that 
it is the line response to a 5-function continuum flash. 

As is the standard practice in the field of line rever- 
beration studies, I identify the time lag between the two 
light cu rves with the peak in their cross-correlation function 
(CCF) dBlandford fc McKee 1982| ). 



Figure 1. The transfer functions ^(i) of various lines as calcu- 
lated for model C. ^ is normalized to 1 at t = 0. 



4 RESULTS 

The present calculations demonstrate the general behaviour 
of the BS system by describing in detail a particular model 
that was investigated in paper II. This model (model C) is 
characterized by a relatively low Lion/A/bh ratio (3.6x10** 
erg s"'^ / 1.9x10*Mq), a BS fraction in the stellar popula- 
tion that decreases like r^^ and a relatively low hydrogen 
density at the base of the wind (3.8x10** cm~^). I calcu- 
lated the transfer functions of various emission lines, which 
are displayed in Fig. ^, by inserting thephotoionization re- 
sults for model C into equations ^ and ^ 

There is a distinct difference between the transfer func- 
tions of the high ionization lines He ll A1640 and O vi A1035, 
and that of the low ionization line Mg ll A2798. This refiects 
the fact that in this model, the high ionization lines are 
formed at the inner BSs, where the ionizing fiux is high, 
whereas the low ionization lines are formed further away 
from the continuum source, where the lowered ionizing flux 
allows the existence of extended, partially ionized regions in 
the BS envelope. This trend resembles what is found in line 
reverberation campaigns, namely that the high ionization 
metal lines have shorter time lags than the low ionization 
metal lines. 

The time lag of a line increases with the width of its 
transfer function. However, the relative widths of the trans- 
fer function give only a rough indication of the rank order 
of the line lags. The time lag of an extended BLR geometry, 
such as that of model C, depends also very strongly on the 
time-scale of the ionizing continuum, because the line CCF 
is the convolution of 9 with the continuum auto-correlation 



function (Pension 1991; Netzer 1991). I check this here by 



calculating the response of the lines to Gaussian shaped con- 
tinuum flashes of various durations (parameterized by their 
standard deviation a). Figure |^ shows two such flashes with 
duration of a = 10 and 100 d, respectively, and the re- 
sponses of the Lya and MgllA2798 lines. As expected, the 
narrow, ct = 10 d flash induces a line light curve very similar 
in shape to its transfer functions and consequently, the two 
light curves respond quickly at lags of 20 d and 32 d, respec- 
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Figure 2. The response of the model C Lya (bold line) and 
MgnA2798 (dashed bold line) lines to narrow (tr = 10 d) and 
wide (cr = 100 d) Gaussian flashes (shaded areas), which are 
centered on t = 0. 



lively. On the other hand, the wide, o — 100 d flash results 
in noticeably larger time lags of 84 d and 120 d, respectively, 
and a larger lag difference between the two lines. 

The response times of various lines as a function of the 
Gaussian pulse width are summarized in Fig. ^. Note that 
even for pulses as wide as ~ 1 year, most of the line lags 
continue to increase with the continuum time-scale. Only 
the responses of H/3 and Ovi A1035 appear to begin to sat- 
urate. The different rates of time lag increase result in the 
surprising effect that not only the time-lags, but also their 
rank order, changes with increasing continuum variability 
time scale. For example, the H/3 line, which lags after the 
OviA1035 and HellA1640 lines when the ionizing contin- 
uum varies quickly, precedes these lines when the contin- 
uum varies slowly. In some lines the time lag is a steep func- 
tion of the continuum variability timescale. For example, the 
HellA1640 Une time lag increases by a factor of 6 when the 
width of the continuum flash is increased by a factor of 18. 
Although I demonstrated these trends only in the case of 
one particular BS model and single Gaussian flashes, I ex- 
pect these results to hold also in the general case. This is 
because all BS models have extended BLR geometries and 
all realistic continuum light curves can be approximated, to 
some degree, by a superposition of Gaussian flashes. 

A realistic ionizing continuum is made up of many peaks 
of different widths. I model this by using the recently ob- 
tained long light curve of PG0844 (Kaspi 1996, private com- 
munication). Figure ^ shows the smoothed observed PG0844 
optical continuum and the response of the HellA1640 and 
MgllA2798 lines. The lags were calculated without the flrst 
100 days of the line light curves, which is roughly the time 
it takes ^' to fall to half its initial maximal value. At times 
shorter than that, the line light curve depends strongly on 
the preceding, unobserved continuum and cannot be calcu- 
lated reliably. Note the difference between the lines response 
to the wide second continuum peak and their response to the 
narrow third peak. This is the pattern of behaviour that was 
suggested by the simulations with single Gaussian flashes 



Figure 3. The response of various model C lines to Gaussian 
shaped continuum flashes, as a function of the Gaussian standard 
deviation cr. 



Line 


lag (days) 


rank order 


hya 


70 


3 


Ha 


75 


4 


H/3 


60 


1 


He I A5876 


70 


3 


He 11 A1640 


65 


2 


C ni A977 


75 


4 


C ni] A1909 


80 


5 


C IV A1549 


80 


5 


NvA1240 


70 


3 


OviA1035 


60 


1 


Mg II A2798 


95 


6 



Table 1. The calculated lags (CCF peaks) of Model C for the 
smoothed continuum light curve of PG0844, shown in Fig. ^ The 
CCF was calculated without the flrst 100 days of the line light 
curve. The CCF is calculated with a 5 day resolution. The rank 
order of the lags (1 is the shortest) is also displayed for conve- 
nience. 



and b y the simulations of Perez, Robinson & de la Fuente 
Such a behaviour w as actually observed in the 1989 



( [1992| ). 



NGC 5548 lUE campaign ( |Netzer fc Maoz 1990| ). Table § 
lists the lags of all the lines calculated here. As with the 
Gaussian flashes, there is an anti-correlation between the 
time lag and the degree of ionization. 



5 DISCUSSION 

In papers I and II, it was shown that the BS model can 
successfully account for two main features of AGN spectra: 
the emission line ratios and their profiles. Here, I study the 
line reverberation of the BS model and check whether it can 
also account for this aspect of AGN properties. I focus on 
the two trends that emerge from the available observational 
results: the positive correlation between the AGN luminos- 
ity and the typical time lag scale and the anti-correlation 
between the line time lag and its degree of ionization. 
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Figure 4. The response of model C He ii A1640 and Mg ii A2798 
lines to the smoothed 4759A (rest frame) continuum of PG0844 
(Kaspi 1996, private communication). In order to make the com- 
parison between the light curves easier, the displayed line light 
curves are shifted to the continuum mean and their variation is 
scaled up to the continuum r.m.s variation. 



Several difficulties stand in the way of a direct compar- 
ison of the BS model predictions with observed line rever- 
beration results. From the observational side, the picture is 
still unclear since the available sample of adequately sampled 
AGN is small and spans only a limited range in luminosity. 
A reliable statistical analysis of the time lags is notoriously 
difficult, and this is compounded by the fact that the quoted 
time lags in the various campaigns were derived by a vari- 
ety of statistical procedures, each with its own biases. The 
strong dependence of the time-lag on the characteristic con- 
tinuum variability is another factor that has yet to be taken 
into account. 

From the theoretical side, there is the difficulty that 
there is some degeneracy among the BS model parameters 
that makes it possible to have similar emissivity maps for 
different luminosities (see section |^ , and therefore the lumi- 
nosity of this AGN model is not well defined. This was not a 
serious problem for comparing the BS model line ratios and 
profiles with the the observed ones, because observationally, 
these do not appear to b e strongly correlated with the AGN 
luminosity ( Netzer 1990 ). This is not the case with the time 
lag results, as there is an ind ication of a correlation betwee n 
the time lag and luminosity (Maoz 1994; Kaspi et al. 1996b). 

To date, published time lag results, which are base d on 



al. 1995D ; NG C 4151 ([Kaspi et al. 19964 [Maoz et 



); NGC 322 7 (Salamanca et al. 1994); NGC 3783 



( gtirpe et al. 1994a ); AKN 120 ([Peterson et al. 1988|; [Peter- 
son fc %iskeU 1991); MRK 279 ([Maoz et al. 199C|; |Stirp~ 



the Ba lmer hues , exist only f or 9 AGN (NGC 5548 ([Ko- 
rista et 
al. 199 



al. 199^ b|); NGC 3 516 ([Wanders et al. 1993[); MRK 590 (^ 
terson ( t al. 1993| ); NGC 4593 ( [Dietrich et al. 1994| )). These 
AGN have estimated luminosities^ in the range 6.5x10''^ to 
1.8x10''* erg s~^ and H/3 time lags in the range 4 to 28 d. 



This sample suggests a weak positive correlation between 
the luminosity and the H/3 time lag. A major improvement 
in testing this correlation is offered by new results from a line 
reverberation campaign on a sample of QSO with luminosi- 
ties from 4 x 10** toSxlO*** erg s ^ . Preliminary estimates of 
the time lags confirm that the correlation is indeed real and 
yield H/3 time lags of approxiamtely 100 d for the luminosity 



range 2 to 5x10*^^ erg s"^ (iKaspi et al. 1996b|) 



Because it is possible that some of the correlation be- 
tween the time lag and luminosity is caused by a correlation 
between the continuum variability time scales and the lu- 
minosity, I use in the simulations the optical continuum of 
PG0844 (one of the QSO in this new sample), which is close 
in its luminosity to that of model C. In doing so, I am as- 
suming that the ionizing light curve is similar to the optical 
one. The luminosity of PG0844, when estimated by the same 
procedure as above, is 4.7x10** erg s~^ and that of model 
C is 2.3x10** erg s~^. I thus expect the lag of model C to 
be intermediate between those of the observed Syl sample 
(lags 30 d for L ^ 2x10** erg s~^) and those of the new 
QSO sample (lags ~ 100 d for L ~ 4x10*^^ erg s"^). The 
model H/3 and Ha lags of 60 and 75 d, respectively, (Ta- 
ble |l|) indeed place it in this intermediate range (as do the 
time lags of all the other lines). Thus, the time lag of BS 
model C is consistent with these observations. 

The BS model time lags display to some extent the 
trend for an inverse correlation between the lag and the 
degree of ionization ('time lag line segregation'), which 
is observed in a few low- luminosity Syls. Two groups of 
lines can be identified among the 11 broad emission lines 
which were calculated. The relatively high ionization species, 
HellA1640, Nv A1240 and Ovi A1035, have narrower trans- 
fer functions and shorter lags than the relatively low ioniza- 
tion species, CiiiA977, Ciii]A1909 and MgiiA2798 (Fig. 
and table ^). There are two notable differences between the 
model results and the observations. The C iv A1549 has a 
large time lag, which is comparable to those of the low ion- 
ization lines and the H/3 has a short time lag, comparable 
to that of the high ionization lines. 

There are reasons to suspect that line reverberation 
time lags of low-luminosity AGN can not be simply ex- 
trapolated to higher luminosities. The results presented in 
Fig. ^ suggest that any detailed ordering of the line time lags, 
especially those with short time lags, strongly depends on 
the typical variability time scales of the ionizing continuum. 
This may explain the seeming ClvA1549 discrepancy. The 
H/3 discrepancy may also be related to the "Lya/H/3 prob- 
lem". All photoionization models of the BLR consistently 
under-predict the Balmer lines emission. This is one of the 
outstanding problems in AGN research, which remains un- 
solved despite extensive efforts. The most likely explanation 
is a combination of reddening and inaccurate treatment of 
the line transfer, which arises from the local nature of the 
escape probability method. This simplified approach may be 
inadequate for the extreme Lya and Ha optical depths that 
are typical of the BLR. Detailed line transfer calculations, 
which deal separately with the different frequency bins in the 
line profiles, are needed (see Netzer 1990 for a review). Such 



t The luminosity is integrated in the range 0.1-1/i, based on the 



measured flux at 
continuum F,, oc u' 



5000 A and the assumption of a power-law 
1/2, Ho = 75 km s-i Mpc'^ and go = 1/2. 
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calculations are beyond the scope of this paper. This situa- 
tion casts doubt on the reliability of the calculated model H/3 
time lags. Note, however, that the agreement between the 
model time lag scale and the observed time lag-luminosity 
correlation does not hang only on the Ha and H/3 time lags 
but is also supported by the time lags of the other lines. 

The time lag line segregation of the BS model is not 
as pronounced as that found in the NGC 5548 campaig n 



(Clavelet al. 1991; Peterson et al. 1991; Dietrich et al. 1993) 



A similar trend is also seen with respect to the line profiles. 
The observed profiles display line to line differences ('profile 
line segregation'), which is more pronounced than that of 
the model. This may indicate that the BS fraction is even 
more centrally concentrated than the assumed /bs oc . 

This rough agreement of the model with the observa- 
tions reinforces two of the conclusions that followed from 
the line profile study in paper H. The line segregation in the 
profile widths and time lags is yet another indication that 
the BSs can not be Comptonized (paper I). Comptonization 
implies that the BS wind boundary is determined by the 
ionizing flux and is at a constant ionization parameter, ir- 
respective of the distance from the black hole. This leads to 
similar emissivity maps for all lines and hence to similar pro- 
files and lags, contrary to what is observed. The line profiles 
and lags also require that the fraction of BSs in the stellar 
population, /bs, fall off at least as . A less concentrated 
BS distribution results in significantly larger lags. For ex- 
ample, BS model B of paper II, which is identical to model 
C apart for having /bs = const., yields a H/3 lag of 210 d for 
the PG0844 continuum. This is in clear contradiction with 
the observations. 

In this paper I studied one particular BS model (model 
C), which is specified by a particular choice of the AGN 
dynamical age, the BS envelope structure and the ionizing 
spectrum. Other choices of these parameters will result in 
other model predictions. This work is part of a series of 
feasibility studies aimed at checking whether one BS model 
can simultaneously reproduce the main broad line properties 
of a typical AGN. Model C was shown to reproduce the line 
ratios and line profile, and was therefore also in the focus of 
this line reverberation study. Direct comparison of model C 
with observed line reverberation time lags must await future 
monitoring campaigns of bright Syls. The question whether 
a BS model of a low-luminosity Syl is consistent with the 
detailed line reverberation data available on such AGN is 
part of the larger question of whether the BS model can be 
extended to the full range of observed AGN luminosities. 
This will be explored in a future paper. In the meanwhile, 
I conclude that BS models, similar to those discussed here, 
are consistent with current line reverberation data on AGN 
in the bright Syl / low-luminosity QSO range. 



REFERENCES 

Alexander T., Netzer H., 1994, MNRAS, 270, 781 (Paper I) 

Alexander T., Netzer H., 1996, MNRAS, submitted (Paper II) 

Bailyn C. D., 1995, ARA&A, 33, 133 

Blandford R. D., McKee C. F., 1982, ApJ, 209, 214 

Clavel J. et al., 1991, ApJ, 366, 64 

Dietrich, M. et al., 1993, ApJ, 408, 416 

Dietrich, M. et al., 1994, A&A, 284, 33 

Edwards A. C, 1980, MNRAS, 190, 757 

Kaspi S. et al., 1996a, ApJ, 470, 336 

Kaspi S., Smith P. S., Maoz D., Netzer H., Januzzi B. T., 1996b, 

ApJL, 471, 75 
Kazanas D., 1989, ApJ, 347, 74 
Korista K. T., et al., 1995, ApJS, 97, 285 

Krabbe A., Genzel R., Drapatz S., Rotaciuc V., 1991, ApJL 382, 
L19 

Maoz D. et al., 1990, ApJ, 351, 75 
Maoz D. et al., 1991, ApJ, 367, 493 

Maoz D., 1994, in Gondhalekar P. M., Horne K., Peterson B. 
M., eds, Reverberation Mapping of the Broad-Line Region in 
Active Galactic Nuclei. ASP, San Francisco, pp. 95-110 

Mathews W. G., 1983, ApJ, 272, 390 

Murphy B. W., Cohn H. N., Durisen R. H., 1991, ApJ, 370, 60 
(MCD) 

Najarro, F., Kudritzki, R. P., Krabbe, A., Genzel, R., Lutz, D., 

Hillier, D. J., in Gredel R., ed., The Galactic Center, in press 
Netzer H., 1990, in Blandford R. D., Netzer H., Woltjer L., Active 

Galactic Nuclei. Springer- Verlag, Berlin, pp. 67—134 
Netzer H., 1991, in Duschl W. J., Wagner S. J., Camerzind M., 

eds. Variability of Active Galaxies. Springer- Verlag, Berlin, 

pp. 107-116 
Netzer H., Maoz D., 1990, ApJL, 365, L5 
Penston M. V., 1988, MNRAS, 233, 601 

Penston M. V., 1991, in Miller, H. R., Wiita P. J., eds. Variability 
of Active Galactic Nuclei. Cambridge University Press, Cam- 
bridge, p. 343 

Perez E., Robinson A., de la Fuente L., 1992, MNRAS, 255, 502 

Peterson B. M. et al., 1988, PASP, 100, 18 

Peterson B. M. ct al., 1991, ApJ, 368, 119 

Peterson B. M. et al., 1993, ApJ, 402, 469 

Peterson B. M., Gaskell C. M., 1991, ApJ, 368, 152 

Salamanca I. et al., 1994, A&A, 282, 742 

Scoville N., Norman C, 1988, ApJ, 332, 163 

Stirpe G. M. et al., 1994a, ApJ, 425, 609 

Stirpe G. M. et al., 1994b, A&A, 285, 857 

Wanders I. et al., 1993, A&A, 269, 39 

Wanders I. et al., 1995, ApJ, 453, L87 



Acknowledgments 



I am very grateful to Hagai Netzer for his advice and 
careful reading of the draft. Shai Kaspi's help and permis- 
sion to use unpublished results are much appreciated. This 
research was partially supported by the Israel Science Foun- 
dation administered by the Israel Academy of Sciences and 
Humanities and by the Jack Adler chair of Extragalactic 
Astronomy at Tel Aviv University. 



© 0000 RAS, MNRAS 000, 000-000 



